Rashba spin splitting in biased semiconductor quantum wells 
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Rashba spin splitting (RSS) in biased semiconductor quantum wells is investigated theoretically 
based on the eight-band envelope function model. We find that at large wave vectors, RSS is 
both nonmonotonic and anisotropic as a function of in-plane wave vector, in contrast to the widely 
used linear and isotropic model. We derive an analytical expression for RSS, which can correctly 
reproduce such nonmonotonic behavior at large wave vectors. We also investigate numerically the 
dependence of RSS on the various band parameters and find that RSS increases with decreasing 
band gap and subband index, increasing valence band offset, external electric field, and well width. 
Our analytical expression for RSS provides a satisfactory explanation to all these features. 

PACS numbers: 71.70.Ej, 73.21.Fg 
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Rashba spin splitting (RSS) due to structure inver- 
sion asymmetry has attracted intense interest recently, 1 
due to its potential application in spintronic devices^*^ 
The origin of RSS, however, has been a controversial is- 
sue for over 30 years. 5 ! 6 ! 7 ! 8 ! 9 ! 10 . Recently, it was shown 
that RSS in the conduction band (CB) was determined 
by the total electric field in the valence band (VB)»ii 
Approximate analytical expressions for RSS were also 
derivecU2ii^. These analytical expressions are helpful for 
understanding RSS. However, they only give rough esti- 
mates of the magnitude of RSS. Most of the features of 
the numerical results for RSS, 12 ! 1 ?' 14 ! 15 ' 16 ! 17 ! 18 e.g., the 
nonmonotonic behaviors at large wave vectors^ and the 
dependence on the subband and well widthi^ can not 
be explained through these expressions. Further, in or- 
der to control the magnitude of RSS, it is important to 
understand the dependence of RSS on the various band 
parameters. However, to the best of our knowledge, this 
issue has not been clearly discussed. 

In this Rapid Communication, we present a theoretical 
study of RSS in biased quantum wells in the framework 
of eight-band envelope function model. The full numeri- 
cal results show that at large wave vectors, RSS is both 
nonmonotonic and anisotropic as a function of in-plane 
wave vector ky , in contrast to the widely used linear and 
monotonic model. We derive an analytical expression for 
RSS, which can correctly reproduce such nonmonotonic 
behavior at large wave vectors. Our analytical expression 
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shows that such behavior comes from the energy disper- 
sion of the subband. We also investigate numerically the 
dependence of RSS on the various band parameters and 
find that RSS increases with decreasing band gap and 
subband index, increasing valence band offset (VBO), 
external electric field, and well width. Our analytical 
expression gives better descriptions to all these depen- 
dences than the previous analytical expressions i 2 * 1 ^ 
In the basis 
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for [001] orientated heterostructures, the eight-band 
Hamiltonian H = i/k=o + + V, where V = 
eFz is the external electric field-induced potential, 
H k=0 =diag{£; c , E C ,E V ,E V , E v , E V ,E V - A, E v - A} is 
the band-edge Hamiltonian, and the upper triangle of 
Hk reads 
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Here 71, 72, 73, k are modified Luttinger parameters, 
(7c — 1) describes the remote band contribution to the CB 
effective mass, k z = —id/dz, and k x ,k y are c-numbers. 

Neglecting the off-diagonal elements in the valence 
bands (VB's) and eliminating the VB components of the 
envelope function, the effective CB Hamiltonian is ob- 
tained asi£ 

h 2 

ffeff(k||) = E c (z) + V(z)+h ^—r k + a (z)(k ]] xe z )-a ) 
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where E p = 2moP^/fi 2 and m*(z) is the effective mass 
given by 
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with U ih (z) = E-Hxk, U so (z) = E-H so , H lh = E v + V + 
P~Q, and H so = E v — A + V + P. Here E is the eigenen- 
ergy and the operators k z in P and Q have been replaced 
by ir/W (W is the well width). The Rashba spin-orbit 
interaction strength a$(z) = h 2 / (fSmo)d^f(z) / dz, where 
j(z) = E p [l/Uih(z) — l/U so (z)]. It is interesting to no- 
tice that RSS at small wave vectors comes from the cou- 
pling to the light hole and spin-orbit split-off VB's, while 
the heavy hole bands do not contribute. This is because 
the basis functions for the electron S |, S J, and those 
for the heavy hole |3/2,±3/2) are spin eigenstates, while 
the k • p interaction between the CB and VB's is in- 
dependent of spin. The last term of Eq. © leads to 
spin-dependent boundary conditions. To obtain an an- 
alytical expression for RSS, we neglect its influence on 
the envelope functioni 2 ^ but we keep it in the effec- 
tive Hamiltonian (we have checked numerically that this 
approximation would not change the qualitative behav- 
ior of the resulting analytical expression at large wave 
vectors), then the RSS of the n-th subband is given by 
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AE^(k {l ) = ^-E p eFk {l J dz \F n (z)\ 2 (U u 2 -U- 2 ). 

(4) 



Here F n (z) is the envelope function of the n-th subband 
along the z axis, zp — Zj ± + , and {zj} denote the z 

coordinates of the interfaces. AEn\li\\) can be viewed 

as the Tg and T7 VBO's-induced interface electric field 

(2) 

contribution, while AEV\ (k||) is roughly proportional to 
the external electric field. For small V(z) compared with 
the band gap, we obtain A£"„ (k|| ) analytically as 
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(If) Here Pi — j, eT .dz\F n (z)\ is the probability of the 
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electron in the j-th layer, w lh 
dropped] are for the j-th layer. Further, if 7 C and the 
wave function penetration into the barriers are neglected, 
and U^, Ui Q are replaced by l/E g and l/(E g + A) (E g 
and A refer to the well material), respectively, then we 
recover the result of Ref. 12, 
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where m* is the CB effective mass of the well. We notice, 
however, that a factor of 3/2 is missing in the definition 
of m* in Ref. 12. Further, Eq. © is invalid for narrow- 
gap semiconductors or narrow quantum wells, where the 
subband energy is comparable to the band gap. 

From the above discussions, we see that RSS comes 
from (i) spin-dependent kinetic and potential energy; (ii) 
expectation value of the total electric field (including the 
external and interface electric fields) in the Ts and Tj 
VB's; (iii) variation of band parameters across the inter- 
faces. The above analytical expressions show that RSS 
is determined by the total electric field in the Tg and 
VB's, in agreement with Ref. 11. As a result, we see that 
the Ando argument^ fails due to an incorrect assump- 
tion that RSS is proportional to the total electric field 
in the CB. This deepens the previous argument*^ that 
the failure of the Ando argument is caused by the spin- 
dependent boundary conditions or the vanishing barrier 
penetration. 

To estimate the validity of our analytical expressions, 
we solve the 8x8 Hamiltonian numerically for a biased 
CdTe/Hgi_ x Cd x Te/CdTe quantum well and compare 
the full numerical solutions with the analytical results 
in Fig. 1. The band parameters can be found in Refs. 
14 and 20. First, we see that our analytical expressions 
[Eqs. Q and 10] agree better with the numerical re- 
sults than the previous analytical expression^ [Eq. J()J)] 
does at small kn. Second, RSS begins to decrease for 
fen larger than a critical value fco. This nonmonotonic 
behavior is correctly reproduced by our analytical ex- 
pression, while the previous analytical expression only 
gives a linear behavior. The decrease of RSS at large 
fc|| arises as follows. With increasing fc|| , the subband en- 
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and becomes comparable to ^(HgCdTe) when fc|| w fco 
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FIG. 1: RSS of the lowest subband vs. fcy for £ g (HgCdTe) = 
(a) 0.8 and (b) 0.2 eV. Solid lines denote the full numerical 
results. Dashed (short-dashed) lines are obtained from Eqs. 
© and © [Eq. ©]. The two contributions AE {1) and AE (2) 
are denoted by dash-dotted lines. We take the external elec- 
tric field _F=100 KV/cm, well width W=20 nm, VBO=0.3 
AE g , where AE g =E g (CdTe) - £ s (HgCdTe). 



The further increase of fc|| leads to the decrease of k^(z) 
and k\\(U^ 2 - C/ s ~ 2 ) in Eqs. © and ©. Consequently, 
RSS begins to decrease when fen > ko. Further, when 
i? g (HgCdTe) decreases from 0.8 to 0.2 eV, the critical 
wave vector kg decreases, since a smaller ko is required 
for E n (k Q ) to become comparable to i? 9 (HgCdTe). Here 
we see that the nonmonotonic behavior of RSS at large fey 
comes from the energy dispersion of the subband. Phys- 
ically, the coupling to the VB's and, consequently, RSS 
are reduced by the increasing energy difference between 
the CB and VB's. 

Third, RSS is dominated by the mean external electric 
field contribution AE^ 2 \ This trend becomes more pro- 
nounced when E g (HgCdTe) is decreased from 0.8 to 0.2 
eV. It can be understood since the interface contribution 
AE^ is roughly proportional to 1/E g through l{zf) in 
Eq. ©, while the mean external electric field contribu- 
tion AE^ is roughly proportional to 1/E 2 through U^ 2 
and U~ 2 in Eq. ©. Fourth, both AE^ and AE^ in- 
crease significantly when E g (HgCdTe) is decreased from 
0.8 to 0.2 eV, in agreement with Eq. ©. Physically, 
the increasing RSS is caused by the increasing coupling 
between the CB and VB's with decreasing band gap. 

In the above, we see that our analytical expression 
gives a better description for RSS than the previous ana- 
lytical expression. Next, by numerically solving the 8x8 
Hamiltonian, we investigate the dependence of RSS on 
the band gap E g (HgCdTe), VBO, subband index, exter- 
nal electric field, and well width. The previous analytical 
expressions^ 2 *^ can only explain the dependence of RSS 
on the band gap and external electric field, while our an- 
alytical expression can explain all the dependences, as we 
shall show below. 

In Fig. 2, we plot the RSS of the lowest three sub- 
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FIG. 2: RSS of the lowest three subbands (the subband in- 
dices are indicated in the figure) vs. ku along [100] (solid 
lines) and [110] (dashed lines) for £ 9 (HgCdTe)= (a) 0.5 and 
(b) 0.1 eV and different VBO's, with fixed well width W=20 
nm and external electric field -F=100 KV/cm. 




FIG. 3: RSS for the lowest subband vs. k x and k y , with 
^=100 KV/cm, £ 9 (HgCdTe)=0.5 eV, and VBO=0.3 AE g . 



bands for different band gaps and VBO's. First, in ad- 
dition to the increase of RSS with decreasing band gap, 
it also increases with increasing VBO. This can be un- 
derstood from Eqs. © and JSJ, because the subband 
energy E n (k\\) decreases with increasing VBO. As a re- 
sult, "f{zf), U^ 2 , and U~ Q 2 in Eqs. Q and (0 increases. 
Second, RSS decreases with increasing subband index, 
due to the increase of the subband energy i5„(k||) and 
the decrease of the asymmetry of the envelope function 
at the two interfaces, because the orthogonality require- 
ment between different cigenstates serves as an effective 
repulsive force, which reduces the potential asymmetry 
produced by the external electric field. 
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FIG. 4: RSS for the lowest subband vs. the electric field 
and well width, with fe||=0.1 nm^ 1 , J3 9 (HgCdTe)=0.5 eV, 
VBO=0.3 AE g . 

Finally, RSS is isotropic at small fen but anisotropic 
at large fc|| . This interesting behavior is in contrast to 
the current understanding that RSS is always isotropic. 
From Fig. 3, it can be seen that RSS has a four-fold 
anisotropy in the k|| space. This anisotropy comes not 
from the macroscopic potential, but from the C4 V sym- 
metry group of the quantum well structure (neglecting 
the bulk inversion asymmetry). 

The dependence of the RSS at fen =0.1 mif 1 on the 
electric field and well width is shown in Fig. 4. The RSS 
increases almost linearly with increasing electric field, in 



agreement with Eq. J5J , while it increases with increasing 
well width and saturates at large well width. The latter 
can be explained through the dependence of "f(z), U^ 2 , 
and Z7~ 2 on the subband energy 25 n (kii). With increas- 
ing well width, E n (kn) decreases due to the decrease of 
the confining energy E n o, such that 7(2), U-^ 2 , U^ 2 and, 
therefore, RSS increase until E n g vanishes and ^(kn) 
approaches a constant value H 2 k 2 /(2m*). Afterwards, 

j(z), U^ 2 , and U~ 2 do not vary appreciably and RSS 
saturates. This behavior is quite different from that of 
asymmetric AlAs/GaAs/Alo.isGao.ssAs quantum wells, 
where RSS shows a peak and then decreases with increas- 
ing well width. 12 It was argued that such behavior comes 
from the competition between confinement and barrier 
penetration. Using our analytical expression, however, 
the origin of such behavior is transparent. That is, in- 
creasing the well width leads to two competing effects: 
the decrease of the subband energy (which increases RSS) 
and the asymmetry of the envelope function at the two 
interfaces (which decreases RSS). 

In summary, based on the full numerical solutions to 
the eight-band envelope function model, we have found 
that at large wave vectors, RSS is both nonmonotonic 
and anisotropic, in contrast to the widely used linear and 
isotropic model. We have derived an analytical expres- 
sion, which can correctly reproduce such nonmonotonic 
behavior at large wave vectors. It shows that the non- 
monotonic behavior comes from the energy dispersion of 
the subband. We have also investigated numerically the 
dependence of RSS on the various band parameters and 
found that RSS increases with decreasing band gap and 
subband index, increasing VBO, external electric field, 
and well width. Our analytical expression gives better 
descriptions to all these dependences than the previous 
analytical expressions. 
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